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Abstract — We build upon and derive a precise far from carrier 
phase noise model for radiation pressure driven opto-mechanical 
oscillators and show that calculations based on our model 
accurately match published phase noise data for such oscillators. 
Furthermore, we derive insights based on the equations presented 
and calculate phase noise for an array of coupled disk resonators, 
showing that it is possible to achieve phase noise as low as -80 
dBc/Hz at 1 kHz offset for a 54 MHz opto-mechanical oscillator. 



In the subsequent sections we propose and analyse an array 
of mechanically coupled disks that is predicted to have better 
phase noise based on insights derived from our model. 

II. Radiation pressure driven opto-mechanical 

OSCILLATORS 

Figure [2] shows a schematic of the device presented in 0. 



I. Introduction 

Opto-mechanical oscillators (OMOs) are unique as they 
feature optical input and output for both power and signal. 
Unlike traditional temperature compensated crystal oscillators 
(TCXOs), their performance is not limited by the / x Q 
product of the material and they are expected to demonstrate 
narrower linewidth and better phase noise. Based on origin 
of self-sustaining oscillations, OMOs can be classified into 3 
main categories: radiation pressure! 1 1 [2 1 [ 3 1, gradient force|4] 
and stimulated Brillouin scattering! 5 1, as shown in Figure [T] 
below. 
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Fig. 1. Illustration of various kinds of opto-mechanical oscillators (6) 

While experimental results have confirmed narrow 
linewidths for OMOs, experimental and theoretical studies of 
phase noise in literature are still limited. In this work, we 
build upon the existing theoretical models and derive a more 
precise phase noise model for radiation pressure driven OMOs 
demonstrated in |2|. We start with an overview of the radiation 
pressure driven opto-mechanical oscillator. Then we present 
analysis of the linewidth and phase noise in these oscillators. 
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Fig. 2. Illustration of setup for radiation pressure driven OMO 

The structure consists of a microtoroid resonator in close 
proximity with a fiber tapered waveguide. Light is fed into 
the waveguide from the tunable laser and the photodetector 
reads off the output. Light couples into the resonator from 
the waveguide through evanescent coupling. The geometry 
of the disk is chosen such that the optical path inside the 
resonator is an integer multiple of the wavelength of the 
light. Due to the high quality factor, the intensity of light 
builds up to very large values enabling strong optomechanical 
interactions. The resulting radiation pressure of the circulating 
optical power exerts a radial force on the microtoroid that 
causes the structure to expand. However, this causes a change 
in the radius and hence the optical path inside the cavity which 
causes a shift in the resonant frequency. This leads to decrease 
in circulating optical power thereby reducing the radiation 
pressure. The restoration of the mechanical deformation sets 
up oscillations in the radial displacement. Usually evanescent 
coupling through a fiber taper is used to couple light into and 
out of the resonator. 

III. Linewidth and Phase noise in OMOs 

The classical Leeson model for electronic oscillators re- 
lates the phase noise of an oscillator to its linewidth. We extend 
this well established theory to the case of opto-mechanical 
oscillators. 



A. Oscillation linewidth 

Mani et al. have characterized the mechanical oscillation 
linewidth as follows[2|: 
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The linewidth is set by fcs (Boltzmann constant), T (ab- 
solute temperature), Q to t (total optical Q), m e ff (effective 
mass of mechanical resonance), Qq (mechanical resonant 
frequency), R (disc radius), T (optical modulation transfer 
function 12), Afirj (intrinsic mechanical oscillation linewidth) 
and M (modulation depth). V represents the finite response 
time taken by the resonator to respond to modulation of the 
optical path length. Intuitively, it should seem that T is a 
low pass filter function of oscillation frequency, i.e. at higher 
mechanical oscillation frequencies its value is smaller. T and 
M are experimental parameters. However, for reasons we 
discuss later, we will write the expression for the linewidth in 
terms of more fundamental parameters describing the system: 
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where r is the amplitude of radial displacement. 

B. Phase noise 

Phase noise is frequency domain representation of random 
fluctuations in the phase of a waveform. An ideal oscillator 
generates a pure sine wave that can be represented as a pair 
of delta functions at the oscillation frequency in the frequency 
domain, i.e. all the signal power is at a single frequency. 
However, due to presence of various noise sources within 
the oscillation loop, real oscillators demonstrate time domain 
instabilities , which appear as random fluctuations in the 
spacing between zero crossings of the waveform from the 
half period of oscillation. This can be interpreted as random 
fluctuations in the frequency of the oscillator, as illustrated in 
Figure [3] 
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where A 2 / 2 is the signal power at the oscillation frequency 
and /a is the linewidth. A/ specifies the offset from os- 
cillation frequency /o. Phase noise, which is specified in 
dBc/Hz (decibels below carrier per hertz), is one metric used 
to quantify the amount of signal power that is contained in 
close offsets from the desired signal frequency. This way of 
quantifying phase noise gives a measure of the spread of 
oscillation power to frequencies offset from the center (carrier) 
frequency Q. 

Phase noise is one of the most important characteristics 
of a self-sustained oscillator. An example of its importance 
can be seen in the use of oscillators as frequency references 
in communication systems: the phase noise of the oscillator 
impacts the bit error rate and security. Also, oscillators do 
not have memory of phase of the waveform and as such are 
unable to restore phase, which leads to accumulation of phase 
deviations caused by phase noise, or drift. Hence it is crucial 
to be able to predict phase noise of an oscillator during the 
design phase. 

We can use Leeson's model to model the phase noise of 
OMOs using the expression of linewidth from equation [2] 
The following equation gives an expression for phase noise 
in both 1/ f 2 and l// 3 regimes. u>i/f3 is the corner frequency 
for l// 3 noise. The l// 3 noise is attributable to slow envi- 
ronmental noise processes Q. 



L(Af) = I0log 10 
= IQlogio 



1 

27 
1 

27 



An 

Ap 

An 



i// 2 



A! 2 ' Af 



I// 3 (4) 



Figure [4] shows that this model gives a close match to the 
measured data. 
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Fig. 3. Illustration of spectral width widening due to noise in oscillators Fi §- 4 - Comparison of model prediction to measured phase noise data |2| 



Thus, in the power spectral density (PSD) of the oscillator, 
the oscillation energy is not located at a single frequency but 
spreads to adjacent frequencies which can be characterized 
by the linewidth of the oscillator. The shape of this spectrum 
for typical oscillators is a Lorentzian distribution, which is 
described by the following equation 



Based on this model, we can draw some insights into 
designing oscillators with low phase noise. As can be seen 
from the equations for linewidth and phase noise, we see 
that if there could be a way to increase the effective mass of 
the system, m e ff, while ensuring that the amplitude of radial 
displacement, r, does not drop so much that — - — ff increases, 



then we can improve the phase noise of the system. However, 
the system becomes stiffer and hence we may need larger input 
power to obtain the same radial displacement. We propose a 
solution to this problem by using multiple disks mechanically 
coupled to each other. Such a configuration allows us to actuate 
multiple disks, thereby maintaining operating power. 

IV. Mechanically coupled opto-mechanical 

OSCILLATORS 

A. Dynamics of radiation pressure driven OMO 

The differential equations governing the dynamics for a ra- 
diation pressure driven OMO have been discussed in previous 
publications CUD. 
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where r is the radial displacement amplitude and Ai n and B 
are input and output optical power intensities as shown Figure 
[5] A is power inside the resonator. 
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Fig. 5. Illustration of opto-mechanical power coupling 

The output power B(t) for a disk can be specified in terms 
of its input power A in {t) and the power inside the resonator 
A(t) as follows 
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where to is the time taken by light to complete one 
circulation of the resonator. When the laser is detuned by Aoj 
relative to the optical mode resonance, the optical detuning 
Aui(t) is given by the following relation 
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B. Mechanically coupled OMOs 

We propose a system of disks mechanically coupled to 
each other as shown in Figure [6] where only one resonator is 
coupled to the waveguide. The figure only shows two coupled 
disks, but one can visualize a linear array of disks coupled in 
a similar fashion. 

The coupling beam length is equal to A mec ^/2, where 
X m ech corresponds to the wavelength of the mechanical oscil- 
lation frequency. This beam ensures strong coupling between 
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Fig. 6. Proposed array of disks 

oscillations of both disks. This configuration increases the 
effective mass for the system while maintaining the oscillation 
frequency. However due to increased stiffness of the system, 
a higher laser power is necessary to maintain the radial 
displacement amplitude. 

To overcome this high power requirement, we can actuate 
multiple disks instead of a single one. Figure [7] illustrates how 
this can be achieved. 
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Fig. 7. Proposed array of disks with multiple disks actuated 

C. Dynamics of mechanically coupled OMOs 

The model for a single disk OMO H) can be extended to 
the mechanically coupled OMO system. The model for a two 
disk system can be interpreted as described in Figure [8] The 
stiffness of the coupling beam is modeled by an equivalent 
spring of spring constant k spring . 
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Fig. 8. Illustration of model for two coupled disk system 

The dynamics of this system can be described by the 
following set of equations: 
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The optical detuning for the two disks can be described as 
follows 
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This system was numerically simulated in MATLAB. The 
half wavelength coupling beam presents infinite stiffness in 
the system [ 10 1. However, since it is difficult to model infinite 
stiffness in time domain numerical simulations, we have used 
a large value for k spr i ng in our simulations. Figure [9] shows 
numerical simulation results for a 2 disk system with only one 



disk coupled to the waveguide. Figure 10 shows numerical 
simulation results for a 2 disk system with both disks coupled 
to the waveguide. As can be seen, the amplitude of radial 
oscillations is higher for the case where both disks are coupled 
to the waveguide. 
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Fig. 10. Radial displacement of both disks - both disks coupled to waveguide: 
amplitude of n = 7*2 = 1 .93pm 



V. Phase noise for mechanically coupled OMO 

Having seen that a system of coupled disks can oscillate 
in synchrony, we can use the phase noise model to estimate 
the phase noise for this system. We can extend the design 
presented in J2 to an array of 100 such disks mechanically 
coupled to each other. Using the model we obtain the following 
data for phase noise at different carrier offset frequencies, as 
shown in the table, showing an overall 20dBc/Hz improve- 
ment due to increased effective mass. 

TABLE II 
OMO Phase noise (dBc/Hz) 
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Fig. 9. Radial displacement of both disks - one disk coupled to waveguide: 
amplitude of r\ = T2 = 1.65pm 

This model can be extended to more than two coupled 
resonators. Table [I] shows radial displacement waveforms for 
a system of five coupled disks after settling. 

TABLE I 

Comparison of stored energy metric m e //r 2 (10 -34 fcg — m 2 ) for 

DIFFERENT NUMBER OF DISKS 



Number of disks 


Stored energy metric 


1 


1.56 


2 


3.72 


3 


4.08 


5 


6.12 



Frequency offset (Hz) 4- 


Single disk(2) 


Proposed design 


100 


-30 


-50 


1000 


-60 


-80 


10000 


-90 


-110 


100000 


-110 


-130 



VI. Conclusion 

We have presented a model for predicting phase noise for 
radiation pressure driven OMOs. Insights derived from this 
model will enable one to design low phase noise oscillators. 
Furthermore, we have proposed and analysed an array of disks 
mechanically coupled to each other and have demonstrated 
better phase noise for the design using the model. The phase 
noise for this array is modeled to be —80dBc/Hz at 1kHz 
offset for a 54M Hz opto-mechanical oscillator. This suggests 
that insight gained from the proposed phase noise model can 
help us achieve improved phase noise performance for chip 
scale opto-mechanical devices. More effort has to be spent on 
understanding the frequency scalability of this model, namely 
understanding how the parameters T and M in equation [T] vary 
with frequency. 
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